SURFACE AND GEODETIC HEASUREBENTS 

(Maryland Univ.) 59 p HC $4.2 



radar optimization for sea 

Final 


CSCL 171 


G3/32 


N75- 18458' 

Unclas 
11790 ; 



TECHNICAL 

RESEARCH 

REPORT 


FINAL TECHNICAL REPORT 


RADAR OPTIMIZATION FOR SEA 
SURFACE AND GEODETIC 
MEASUREMENTS 


National Aeronautics and Space Administration 
Grant No. NGR 21-002-433 



DEPARTMENT OF 
ELECTRICAL ENGINEERING 


UNIVERSITY OF MARYLAND 

COLLEGE PARK, MARYLAND 20742 



6 FINAL TECHNICAL REPORT 

RADAR OPTIMIZATION FOR SEA SURFACE 
AND GEODETIC MEASUREMENTS 


National Aeronautics and Space Administration 


Grant No. NGR 21-002-433, 
• 1/1/74 to 12/31/74 


Principal Investigator 
Dr. Robert O. Harger 


Department of Electrical Engineering 
University of Maryland 
College Park, Maryland 2 0742 



EFFICIENT GEOID AND SEA STATE 
ESTIMATION BY SATELLITE ALTIMETRY 


Robert O. Harger 


August 1974 


Electrical Engineering Department 
University of Maryland 
College Park, Maryland 20742 



Research supported by NASA Grant NGR 21-002-433 



ABSTRACT 


The efficient estimation of geoid and sea state parameters (waveheight 
standard deviation and asymmetry and scattering cross-section density) is 
considered: the optimum processing structures, including maximum 
likelihood estimators, and their accuracy limits are given for a model 
accounting for random surface reflectivity, sea height, aind additive noise, 
and allowing for arbitrary radar system parameters, based on the assumption 
the received signal is a sample function of a normal random process. 

The integral equation associated with the "Gaussian signal in 
Gaussian noise 1 ’ inference problem was solved previously in 11 pulse- 

z' t 

determined" and "threshold" cases and is here solved approximately in 
the quasi-stationary case of pulse resolution small relative to the returned 
signal's extent and exactly in a specific case.. It is shown that the optimum 
processing is generally a mixture of coherent and incoherent integrations 

which may be viewed as a weighted summation of received power of the - 

match- filtered received data. There is a signal-to-noise ratio dependent 
best pulse bandwidth; approximate and exact accuracy limits are given. 

The efficient estimates are generally correlated, the most strong coupling 
between geoid and asymmetry estimates and between wave height standard - 
deviation and reflectivity estimates. A brief exemplification of system 
design usage is given. 



INTRODUCTION 


Radar altimetry has been proposed and is being used to obtain 
measurements, over the sea from a satellite, of geoid and sea state 
parameters which are expected to be highly useful in geodesy, oceanography, 
and meteorology both for scientific uses and applications such as traffic 
routing. The accuracy required of geoid measurement for " dynamic 
oceanography " is about 0.1 meter; since tides, currents, and wind-driven 
waves cause sea surface perturbations much greater than this and since the 
reflection from the sea surface is a random phenomena, the extraction of 
geoid information from the “ smeared" and random returned pulse 
requires (at least implicitly) joint estimation of tides, currents, and 
sea state parameters. It does not appear possible for altimetry itself, 4o 
separate tidal, current, and geoidal variations: exterior data must be 

provided for this. But sea state parameters do directly influence the 
nature of the altimeter’s returned signal from the sea and must be estimated, 
at least implicitly, whether or not desired in themselves, in order to obtain 
a good geoid estimate. - ' 

— — —W-e- -consider here the be st -single -pulse joint measurement of the 
_geoid parameter - the delay of the returned pulse from the mean sea level 
(which here inseparably includes tide and current changes) - and sea state 
parameters which are here taken as the radar scattering cross-section 
density, standard deviation of the wave -height distribution (which is assumed zero 



mean since any nonzero mean cannot be distinguished from the geoid by 
definition) and an asymmetry parameter accounting for the well-known 
asymmetry of ocean waves. Such measurements can then be combined with 
adjunct data - e. g. , laser measurements of orbit - to estimate the geoid; 
there are well-known methods (e.g., least-squares, extended Kalman 
filtering) to do this and they are being so used. Since it appears that 
practical altimeter designs are possible that provide accurate measurements 
even on a single pulse, a linear observation model in a recursive estimation 
scheme may be attractive; a general statement of a recursive estimation 
problem appears to be very involved, hence the emphasis on optimizing 
a single pulse's estimates. 

The present discussion is an extension of [ l] where estimation of geoid 
only was considered. The introductory remarks of [l] with respect to 
theoretical context and desirability of efficient design apply here also and 
we adopt the model for the returned signal and the notation of [l]. In 
particular the signal random process covariance function in the symmetric 
case is .... 

dyq(y) J dx Ph (x)F(t i' T 0 ' y ' _ f ] F( W y "'f 1 (1> _ 

0 -® 

where F is the complex modulation of the transmitted pulse, T q the delay to 

geoid, c the velocity of light, p^ the univariate probability density function 

of wave height which depends on the unknown sea-state parameter a^ f 

q(y) = g (/cR~y) H(/cR Q y) where g Q is the reflectivity density and H is the 

fourth power of the modulus of the far-field antenna pattern, both a function of viewing 

angle from the vertical, and K = 2 tt (k/4irR^) (cR Q /2) where k=2TT/\, X the 



mean wavelength of the radiated pulse and R^ = cT^/2. It has been assumed 
that the reflectivity of the sea is a random field with a reflectivity density 
g(x, y) = g (x, y)/secy(x, y) where y(x, y) is the angle between the normal to 
the gross roughness surface of the sea and the altimeter bores ight direction 
and g^(x, y) is a complex random function of zero mean and spatially 
incoherent - that is, its covariance function is of the form gg( X }» Yj) 

6 (x^-x^, y^-y^). (See Fig. 1 for the geometry of the problem. ) 



For the narrow antenna patterns of interest here, g (/cR y) w g (0) — g » 
generally an unknown parameter* that must be estimated along with the 
unknown geoid delay (T ) and sea state parameter 0^. It is then convenient 
to define q^(y) = H(/cRjy) and rewrite Eq. 1 in the form 

WV = g oJ f(5 - T o ;CT h 2)F<t l -5)F(t z -?) d? (la) 

_ oo 

where the convolution 

f(? -V C h 2) 3 J 0 q l (Tl)p h C f ( S-V T ' )]d11 <2). 

2 2 

Note the unknown parameters (T , a , g ) enter via g f{£ - T ; CT ) - which 

r oho o o n 

is denoted more briefly f^(£). 

As the received, scattered signal is modelled as a sample function of 
a nonstationary, zero mean, complex Gaussian random process - an 
assumption with some theoretical and experimental support - whose 
covariance function depends on geoid and sea state parameters, and is 
received along with additive, white thermal (Gaussian) noise, the likelihood 
ratio - necessary for statistical inference - involves the quadratic functional 
of the received data appropriate to the " Gaussian signal in Gaussian noise 
problem*' [l, 2] having a kernel function determined by an integral equation 
that cannot, unfortunately, be solved in the generality desired here. Two - 
special cases were solved in Tl], the M threshold " case and "pulse 
determined" case, and those results are here extended to the present joint 
estimation problem: however, the pulse determined case would a priori 
seem to be possibly unattractive for sea state estimation. 

*The relation between g and c , the radar scattering cross-section density, 
m2 o ° 

is O = TTg /X . 
o & o 



We show here that in the threshold case the pulse resolution to best 
estimate sea state should be somewhat smaller than to best estimate geoid 
but still of the order of the received (smeared) pulse duration. To adopt 
this result is generally disquieting because there might well .be a best 
“SNR per sample” (resolution element) analogously to the diversity communication 
problem, and if so, a SNR- dependent, possibly finer resolution would be 
dictated. 

Therefore we are led to consider an approximation appropriate to the 
case when the pulse resolution is small relative to the received pulse 
duration as spread by the antenna pattern and wave height distribution. We 
are able to approximately solve the required integral equation in this 
“quasi-stationary " case, observe the structure of the optimum processor, 
and calculate bounds on estimation accuracy - which are attainable asymptot- 
ically by maximum likelihood estimation. As an example application to altimeter 
system design, there is shown to be a best pulse resolution for estimating 
geoid and sea state, dependent upon SNR. These results can be used for 
efficient altimeter system design as exemplified in [l]. 

Also the integral equation is solved exactly under a specific set of 
assumptions; the results generally agree with the results of the quasi- 
_ -Stationary approximations when the signal bandwidth (roughly- see below for 
more precise statement) is relatively large. 

Since the received signal is modeled as a sample function of a random 
..process, the better known theory appropriate to a signal known except for a 
set of parameters (which leads to the matched filter, etc. ) is not applicable here. 
For comparison, the estimate of delay and spread for a non - random signal is 
discussed in the Appendix. 



SUMMARY OF RESULTS 


The data-dependent part of the log-likelihood ratio, which describes 
the operations on the data Z^, te T, performed by the receiver/processor is 
[l,2] the quadratic functional 


: Q < Z)= dt i z(t i> J 

T T 

where the kernel k^jt^) is the solution to the integral equation 


WV + J R . (t l- t 3 )k( W dt 3 = R s (t l-V> VV 6 T - 


(4) 


(3) 


T is the time during which the scattered pulse is received and will be 

assumed large relative to the support of R ; rj is the spectral density of 

s o 

the complex white noise. 

Solutions for k. - We have already solved the integral equation (4) in 
[l] in two special cases. First, in the "threshold case" (see Tl] for 
definition) of, roughly, low signal-to-noise ratio, 

k <vv « V' VVV ' • (4a) 


which gives 


Q(z) - -^q(-t)*p h (- | t)*|F(-t)*z(t)| 


(3a) 


t=T 


which can be calculated by the processor shown in Fig. 2. As discussed in , 

[l]this processor performs both coherent and incoherent integration in an 

intuitively satisfying manner; since the scale of q involves g^ and depends 

on 0 (and the asymmetry parameter introduced below) these matched 
h 

filters will depend on sea state parameter estimates. 



Second, in the 11 pulse determined 11 case (see [l] for definition) where 
the radiated pulse has insufficient bandwidth to resolve the wave height and 
antenna spot smearing, 

k <V t 2 ) = k o F( V T o )F( V T o ) <4b) 

CO 

where k = a/(T| + cfC ), £ the energy of F and or ='K q(y)dy. Now 

r\ n r « 

O 

2 


o F” F 
k 


Q(Z) |F(-t)*Z(t)| 


<3b) 


which, aside from scale, agrees with the results of the threshold case with 
q(_t) = p (- — t) = 6 (t) - reasonably because the transmitted pulse does not 
11 see 11 the waves and antenna spot in this case. The calculation of 3b is 
accomplished by the processor of Fig. 2 with the matched filters of impulse 
responses q(-t) and p^(-ct/2) omitted. Such a processor is therefore 
independent of the sea state parameters (aside from g Q ) which is an advantage 
if only geoid estimation is required. 

We here solve this integral equation in two more cases.’ Thus, thirdly, 
in the u quasi- stationary " case where the extent of F (after matched filtering 
if dispersed modulation is employed by the transmitter) is small relative 
to the distance over which q(t) * t) varies significantly - that is, the gross 

’waves and antenna spot are relatively well resolved - we find that k(t^, t^) 
is~ approximate by a slowly varying impulse response k^, t^-t^). That is, 
for each t^, in a region of the order of the transmitted pulse resolution, k 
looks like the impulse response of a linear, time -invariant filter. Its 
11 quasi- stationary M transfer function is shown below to be approximately 

VV“>= — 2 < 4c > 

^0+ f *M F(UJ) i 



•which is a weighted sum of the received instantaneous power after matched 
filtering to the transmitted pulse. (Here 2rr£^_,/f)^ where is the bandwidth of F 
which, for the specific form (3c), was assumed constant over its spectrum.) 

The calculation of Q(z) is again accomplished by the processor configuration 
of Fig. 3 except the matched filters q and p^ are replaced by a filter of 
impulse response fJt)/[Tl /0+ Mt)] . This form of processing has in fact 
been a starting point for practical designs [ 5] . A sketch of the weighting 
function w(t) sf (t) • r p /$ + f . (t) ] 1 for a crudely approximated f^ is shown 
in Fig. 3 ; exact calculations are easily made for a given f^. 

It should be noted that w(t) is relatively broad in its extent. The weighting 
function does depend on the sea state parameter estimates. 

- . The first case overlaps the second and third cases which are mutually 
exclusive; taken together they lend great weight to using a processor of the 



an envelope detector and filtering (incoherent integration), the exact filter 

shape dependent on signal- to-noise ratio and sea state parameters. Such 

a processor is eminently practical also. Again if only geoid estimation is 

desired this filter may be omitted and satisfactory performance is practically 

achievable T 1 ] without need for sea state parameter estimation. 

Fourthly and finally, if it is assumed that the transmitted pulse 

spectrum is constant over its bandwidth and that q tp is the Fourier transform 

h 

of a rational function, an integral equation equivalent to (4) can be solved 
exactly. The method of solution is well known and has been reduced to 
routine calculation which gets sufficiently involved to impede easy under- 
standing. A relatively simple exponential form for q cp^ was assumed and the 
calculations performed to find k exactly (Eq. 25 below). It can be expected 
in practice that the system parameters will be such (see below) that the 
resulting processor has the realization shown in Fig. 2 where the filters 

q and p " are replaced by a filter transfer function 

n a . 


~ ~ - Bmjj 

k(w) = k e K 1 1 , -«< uj< co . 

o 


(3d) 


a ij 2 

Note that its bandwidth is about P = {r\j2a§ <(?) where ^ is a measure 

of the spread or smear of the returned pulse: this can be seen to be in 

agreement with the result under the quasi- stationary case. 

Maximum likelihood estimates. - By definition the maximum likelihood 

...estimates ("MLE’s") are. those values of the .unknown parameters that 

maximize the likelihood ratio or, equivalently, the natural logarithm of the 

likelihood ratio which here has the form [ - § + Q(z)] ; $ is independent of the 

data and here can be shown [ l] to be independent of geoid parameter T but 

o 

generally will depend on the sea state parameters. * Thus the maximum 


^Finding useful '’closed 1 ' forms for $ corresponds in difficulty to find k, the 
kernel of Q. 



likelihood estimate of T is easily seen to be found by choosing the time of max* 

o 

, imum output of the low pass filter for any and all of the special cases discussed. 
Said filter must be adjusted over a trial set of sea state parameters to 
maximize the output: in principal - but not necessarily practically - this may 
be done "in parallel" on one pulse return. 

Necessary conditions for the MLE’s are found, of course, by setting the 
derivatives of the log-likelihood ratio with respect to each unknown parameter 
equal to zero. We show below that in the quasi-stationary case and when F 
has constant spectrum over its bandwidth 0^,, that these conditions are 


■i 


dt 




C-n /g+f (t)] 

o * 


{ 


PB. 


o # 


(t)]} . 


( 5 ) 


These likelihood equations can have solutions for each stationary point including 
the true maximum and are therefore meaningful in suggesting processors use- 

A 

ful after initial acquisitions to a close approximation of the estimates 0 to 

the true values ^.= (6^ 6.,, 63, 64) = (T o> 0 ^ , g o> 64), & 4 the asymmetry 

parameter. The functions v^t) s (3 fj 3 0 .) T n q /P + fjt)] " Z may be regarded as 

weighting functions and cjre sketched in Fig. 4 . for a crude approximation to 

_a typical f^ when p^ is symmetric. ' 

-Such-likelihood equations have often been used to suggest feedback 

' structures to find estimates by operating on a sequence of received pulses and 

driving certain error signals to zero. For example, the form of v^(t) 

suggests the use of two time " gates", one with the weighting function equal 

to v (t) for -0 <t < o and the other with the weighting function equal to the 

2 



the negative of v (t) for o < t < a+ cj>. Each gate output -weighted 

Cm 

instantaneous power-is summed and the difference formed as an error signal 

which readjusts the position of the gating pair. The integral of v^ft) times 

[n /0+f (t) D is computed separately and applied as a bias to the difference, 
o * 

Similarly for v^ and v^. 

Performance. - It is well known [4] that the covariance matrix of 

efficient estimates - those unbiased estimates of minimum possible error 

variance which are, when they exist, MLE's and whose performance is 

achieved asymptotically by MLE's - is the inverse of the matrix C of 

2 

elements c.. =- E{ 3 l n A(z)/B0 6- } where A(z) is the likelihood ratio 

and f 0.} the unknown parameters to be estimated. Here the {c } can be 
1 

expressed in an integral involving R g and k (eq. 19 below, e. g. ).' For 
example, in the quasi- stationary case we have the approximation. 


J 


ITT ' W. V" a * ■ 

1 1 


= - f 

2rr J 


B fl(a» * fj® )do>. 


he. * 


5 0. 


( 6 ) 


Here ^ = B /(n /0 + S *(©)) and we have assumed | F (U) )| = 2 tt t / 0 = 0 
i? O O x r 

for |tt)| < Qp/ 2 and zero otherwise (B^ = 0^/211 ); also we have made a 

certain approximation with regard to the general shape of - see below. 

Here the Fourier transform of f is 1 ((D) = g K q 1 (tt))cp T u> / (c/2)] . 

exp(-i(JUT ), where cp is the univariate characteristic function (the Fourier 
■ o h 

transform of p, ) of the gross wave height, 
h 

An often used [4] approximation to the gross wave distribution is 
normal; but it is. usually still better approximated by incorporating a 



parameter that accounts for the observed asymmetry. This is conveniently 


done by using a truncated Gram-Charlier series T 3] : ' 

-a 2 § 2 / 2 

V §)=e ■ a-ie 4a ? J ) 

when 0 = E[h )/3 is proportional to the (assumed small) third moment 

of the gross wave height distribution; note 0^ = 0 if p is symmetric. 

We can now compute the { c.^} and invert C to find the desired 


covariance matrix; its diagonal elements are the error variances of the 

h 


x 

efficient estimates (T , a. f g ,0 ){0 A = 8 A /(c/2) , 0 sa^/( c /2)): 

O 1 0 4 4 4a 




Var { a ] = 


and 


1 

1 

l 

"V 

U P 14 2 

M(l) * 

1 

1 

1 

<t 2 
g o 

* . 2 
1 “ P 14 

‘ M(3) ’ 

1 

1 

1 

V 

. 2 
1-P 23 

o Z M(2) 



2 

1 

1 

g o 

■ V 

*V 

M(0) 1 




(7a) 

(7b) 

(7c) 


<7d) 


where 


2 M( 2) 

14 "M(1)M{3) 1 

2 


2 __ M(l) 


23 


M(0)M(2) 


and 


M( k) s ^ Jlu Zk ] ftiD )] 2 dll) , 

The correlation coefficient of the 


(7 e) 
(7 f) 

(7 g) 


7(«) = Kg («))cp (u)/(c/2)}/ g Q . 

1 y 4a 


efficient estimates of geoid delay and asymmetry parameter 0^ is p^; 


the correlation coefficient of the efficient estimates of gross wave height 



standard deviation a- and sea reflectivity density g is p_ ; all other correlations 

O Ct j 

are zero to this first order (in 0^) approximation; also p ^ and are 

independent of 0 to first order.. 

Consider, for example, the error variance of the efficient estimate 

T of the delay T to the geoid. It is uncorrelated with the efficient 
o o 

■■if 

estimate g of reflectivity density though, because additive noise is 

present, its magnitude depends on the true reflectivity density g^; this is 

* 

agreeable. It is also uncorrelated with the efficient estimate o of the 

standard deviation a of the gross wave height distribution - a second 

moment - though its magnitude depends on this measure of (symmetric) 

spread as is agreeable because greater spreading means greater returned 

energy which, if processed properly (as it is), should enhance the geoid 

estimate. Finally it is correlated with the efficient estimate of the asymmetry 

parameter: errors in estimating the asymmetry of the gross wave height 

distribution would be expected to lead to geoid location errors. Finally, we 

2 2 - 

may write M(l) = where 0 f is interpretable as a radius of gyration 

. ** - 2 

measure of the bandwidth of f ( provided the centroid of | f j is zero): 
this too is agreeable - though not with the superficial thought that the 
bandwidth of F, the transmitted modulation, determines geoid measurement _ 
accuracy. 

# * 

The correlation between the efficient estimates a and is also 
agreeable: increased total energy is scattered by increasing either a or 

v e - g - 

2 1 2 2 2 

More specifically, assume |.q(w)l - | q(0)| exp(-<|) w ) so that 



]f(uu)|^= a^exp{- 0 tf^/ 2 v^) where = K|q^(0)| and v - [ 2 (ct +<{> )] . 

2 2 2 -1 
Then one calculates p = 9/15 and = 1/3, leading to (1- P 14 ) = 15/6 

and (1- p ) =3/2, these being the factors by which the various efficient 

estimates* error variances are increased by their being correlated. Or, 

to say it another way, the variance of the geoid delay efficient estimate 

(e. g. ) is increased by 250% by the lack of knowledge of asymmetry. Further, 

after choosing an optimal pulse bandwidth as discussed below, 

Var{T } = 


and 


i£/ 

o 

TVf \ 

'*o*l*F 

' 2 

: f 110 

3TT./2 

Vi e : 

± ( 

T1 °' ) 

Tr/2 ' 

g o“i e F 

b 

/ \ 

tt/T 

^o“i e F 


)i 


A 4 


(8a) 


(8b) 


(8c) 


(8d) 


8o 

Because of the aforementioned approximation to the shape of f + , it 
was of interest to compute the {c^.3 for a specific, reasonable f^ with 
a symmetric gross wave height distinction - still under the quasi- stationary 
approximation and the specific, reasonable assumption of the shape of (Fj. 
The results are given below (eqs. 22), the result of tedious residue 
integrations; for comparison the approximate forms (Eq. b) are also 
evaluated with results shown in Eq. 23. (Both 22 and 23 result after an 
optimum choice of signal bandwidth 0^ - discussed below.) We see that 
agreement is very close. 



Also derived below are the error variances of the efficient estimators 

2 

of geoid delay T and gross wave height spread in the threshold case, 

For example, for a normal gross wave height distribution, |F| constant 


over its spectrum, and |q(u)| 2 - ]q (O)exp(- <f)jU /2)| , then setting V — 


~2 _ 


2 , 2 


8o^ / c + and a = vQ^,, we find 

* , V N 2 1 

Var{T }= ( ~z p — J * 

° VZg o“i e F ; r ll (a) 


Varfa } 

h 

(c/2) 4 


I * 

= ( ~ 

\2e a 


Z V*i e F 


r 22 (£) ‘ 


^ A 

The forms r^a) and r^a) are graphed in Fig. 6 as a function of a. 

* 2* 

The efficient estimates T and a are uncorrelated - as is proved true 

oh 

for any symmetric p^. 

. Finally we recall from [l] the variance of the efficient geoid delay 

$ 

estimate T in the pulse determined case is 
- o 

Var{ T *} = -2-2 ^ • -±r. (10) 

< Zg o“lV A F 


where A is the radius of gyration measure of the bandwidth of F. (When 
F 

|f| is constant over its support then - Q^/2/3 .) 




III. APPLICATIONS 


Modulation bandwidth. - One of the fundamentally interesting questions 
about system design is the choice of transmitted pulse bandwidth to minimize 
the estimation error variances. 

Taking up first the quasi- stationary case where the bandwidth ft ^ of 
the pulse F is such that the *' smeared" returned pulse f is highly resolved, 


we note - Eqs. 7 - that 0 enters only through <£ and each of the error 

F 

variances are inversely proportional to CL : thus ft should be chosen to 

r 

maximize Cl = B„/[t| B / S + g f{0)] where B =ft /2tt. It is easily 
F o F F o F F 

seen that the maximizing B is 

g f L 
o o F 


max 


4g f *n 

o o o 


To interpret this result note that, from Eq. lb, the mean power of the 
received signal at time t } (said mean power will slowly vary in this quasi- 
stationary case) is a s 2 (tj) = R^t^) - , where ^ is the energy of 

the complex modulation (which is twice the energy of the real pulse modulation); 

also the mean power of the complex noise in the bandwidth B is 

2 F °P l 

°n = T1 o B F ' 0Uf result states that B p should be chosen so that the signal- 

2 2 

to- noise ratio G g (0 )/o^ is unity (recall we set = 0 for present purpose). 

Since f(t) generally will decrease as 1 1 | increases from 0 and since we 

- nevertheless can only choose one we may expect that a more exact 

calculation would result in a somewhat smaller B , corresponding to 

4 . 2 2 

a ratio a g (°)/ a n somewhat greater than unity. '■ • r' he -;: 



In fact this last is borne out by examining the results of the exact 
calculation, Eqs. 21e, 21f, 21g, for the symmetric case. Each variance 
depends on a different factor dependent on a = (1 + a/Q^J , only through 
which 0 enters: a can equivalently be chosen, then, to minimize (e. g. ) 

F 

one of these three forms-call them CT; (a), i= 1, 2, 3. It is easily found algebraically 
that a=2 minimizes both Var(a*) and Var(g Q *) and a=(3+/l7)/4 »1. 8 minimizes 

Var(T *); further, a choice of a=2 increases Var(T ) by less than 2%. Thus the 

° 2 2 
choice a= 2 is a very good one; then .2/^ = (2 tt f^CO)^/^)/ 2 ^ 5 -^, = (°)/ a n 

V is seen to be (very nearly simultaneously) optimally 3. 

This choice of a is allowable, provided, roughly, the signal-to-noise 

ratio is large enough so that the quasi- stationary approximation remains 

/v 

valid: for this we must have 1 / 0^ << 9. Since 0/0=3 (for a= 2), 

1/0= 3/a = (3/2)(T^/ae F )0: thus, more precisely, we must have 

T) 

“ 7 “ << 1 . 

F 

As will be seen below, this condition is easily met in practice even for a 
satellite vehicle. 

As will be seen in examples below the optimum choice of a can imply 
very large pulse modulation bandwidths. It is therefore of interest to see 
how the { <t .(a)} depend on a for a greater than the optimum a - that is, 

'bandwidths less than the optimum bandwidth: this dependence is sketeched 
in Fig. 5. To conform to the quasi- stationary approximation only a‘s below 
a bound are permissible: that is, for a given 4tt q* ^p/^g' since < 1, 

only values of (4n a £ F /n o ){p F /cj>) << Hna^/rM are allowable - which is 
to say only a 2 <<(4tt ) 4 1 are allowable. For example, in an application 



discussed below (4TTa£_/r| } is greater than 10 6 so a’s somewhat less than, 

± o . — 

say, 300 are allowable. 

We see that the error variances increase by an order of magnitude as a 

- 1/2 

increases by an order of magnitude: fortunately, roughly, a.— ( 0 ^) so 
that an order of magnitude change in a corresponds to two orders of 
magnitude change in 0^,* 

The best choice of bandwidth 0 F in the threshold case is also examined: 
0 enters the error variances (Eqs. 9) only through y^(a) andy^fa), 

a = vfi , which are graphed in Fig. 6. 

F ...... 

We note that the best choice of Q = a/v (fixed v) is «4/v» degrading 
only a little the accuracy of the geoid estimate and would represent a 
reasonable compromise. • ' • - 



rJ2 2 2 2 2 

As v = 8g /c + » and a, is unknown and to be estimated, we 

h i n 

2 2 2 

may have a problem unless <j> » 8a, /c - that is, unless the antenna 

1 h 

•'spot" determines v. Since, from Fig. 1, too small a bandwidth results in 

serious degradation, we may wish to choose 0 & 4/v . , where v . is 

F mm mm 

the minimum possible v - very likely determined practically by the 

antenna spot (scaled into time and smeared by pointing errors) as the sea 

2 

could be smooth (a^ = 0) occassionally. Presumeably generally better 

performance could be obtained by an adaptive system which utilized an initial 

estimates of a to more optimally set Q on the subsequent pulse, 
n F 

System design. - We briefly indicate how the above results may be. 
used in determining reasonable system parameters. 

Recall (Eq. 22 of [l]) 


or = 


AaQ 

(411) 3 R 2 
o 


(13) 


where A is the antenna aperture area, the scattering cross-section density 
2 

O = TTg /X (a measured and tabulated parameter), and Q depends on 

o o 

aperture illumination function shape* e. g. , Q = 1/6 for uniform illumination. 
Also (App. Ill of I) 

2 

X R - (14) 


<!>! » 


Ac 


Since <j>^ decreases and a increases with increasing physical antenna 
aperture area A, clearly we want A as large as possible - limited by vehicle 
size and stability as discussed in Fl] . 



Suppose certain measurement occuracies are required of geoid and sea 
state estimates. Rescaling to spatial distance instead of time extent by 
setting T = R /(c/2) and c = c^/(c/2) , we have 


and 


Var (R ) = (4. 51) 

O 0 ! & 


* . . . \ 


Sir"* <•*?>*].* 


F h 

2 2 6 

Assume parameters o. = tt , A = IT m , R^ “ 10 m, and Q = 1/6: then 

we calculate a- (3. 84 x 10^) If \ = 2 cm, then cf> c/2 = 25.4m. If 

-21 -7 

T| / 2 = kT = 1. 4 x 10 , then T\ /a- 5. 38 x 10 . Thus (Eqsw 8) 

o o o 

Var (R ) = 1. 56 x 10 3 / 1 

o i? / . 

and 

Var (a*) = 8.03 xlO’V^a 2 
h r n 

2 2 

provided that a <<(25.4) , as would be true even in the roughest possible 

_h. 


seas 


. If */ Var(R ) = 10' 1 m. accuracy is specified then at least 
o 


= 1. 56 x 10 * Joule is required, which is rather large; then Vvar(a^ ) = 

2.27/g m. which of course will exceed any specified value for sufficiently 
h 

small o. « 
h 

2 2 

If R is chosen so that c_ (0)/a = 3, then we find B ** 1 KMHz. , an 

F b n r 

~ imp ractically large value. But in order even to be in the fine resolution, or 
quasi-stationary case, the resolution l/B^, must be small relative to 
= 0.17 usee, implying > >5. 9 MHz. 

Suppose a more practical bandwidth of B^ = 33 MHz, which is a 
thirtieth of the optimum bandwidth: then it is easily calculated that now 



a = 10. (We already saw that a = 2 would be increased roughly by 

opt 

/3Q , near optimum a little less.) From Fig. 5 we see that now CH^(a) = 135 
and - (a) = 148, increases from the minimum values of 18 and 27, resp. , 

u 

corresponding to increases in the error standard deviations by the factors 
/l3f>/18 « 2. 74 and /T4l?/27 «2. 34, resp., which would have to be made up 
by changing other system parameters to meet the same specified accuracy. 
This calculation points out that the error standard deviations are relatively 
mildly dependent on pulse modulation bandwidth, within the quasi- stationary 


approximation. 



It is of considerable interest to compare this : 

* 2 2 

achievable performance, Var(T^ ) = (4- 51) (rj a£p)( a + <j>j ) with that 

achievable in the " pulse-determined" opposite extreme pulse bandwidth 

Cl where q is not resolved at all, namely (ft], Eq. 24. b) Var(T Q ) - 
F 

? -2 
6(T] / at ) C2 where (see [l], Eqs. 9 and 10) Q_ is larger than 

(O 2 + . Thus, set O^," 2 = 6 2 (a 2 + c^ 2 ): a general conclusion of Tl] 

was that, in the threshold case, 6 = 4/3 was optimum, resulting in a pulse 

that did not quite resolve the returned signal; adopting this 6 , we see that 

in the pulse determined case the achievable accuracy is at best approximately 

(24) (p fat ){a 2 + <j> 2 ). 

.0x1 

Thus the best performance, measured by error standard deviation, 
in the fine resolution extreme is only about / 24/ 4. $1 ~ 2 times better than 
that obtainable in the pulse determined case. Since the processing required 
in the latter case is markedly simpler and a priori determined, whenever 
only geoid estimation is desired the pulse determined case seems 
markedly preferable: such design examples were considered at length in 

Cl]. 



IV. DERIVATIONS 


Quasi- stationary case. - The data-dependent part of the log-likelihood rs 
which describes the operations on the data Z , t G T, performed by the receiver/ 

v 

processor is [1,2] the quadratic functional 


Q (z)=iJ dtjZttj)/ At 2 z{t z )Mt v t z ) 

■oT T 


( 15 ) 


where the kernel k(t 1# t 2 ) is the solution to the integral equation 



t 3 )k<Wdt 3 R s (t l’V' 


w 


T. 


( 16 ) 


T is the time during which the scattered pulse is received and will be 

assumed large relative to the support of R ; Tj is the spectral density of 

s o 

the complex white noise. 



Suppose- the support of F is small relative to the distance over which 
f changes significantly: then, from (la), 

WV" g o f(t i- V°h 2) * (t rV . ' ' < 17 > 

= f *<V* ( W 

where 9 = F* F , * denoting convolution. 

Note that if in fact f (t) were independent of t then the signal process 
would be wide sense stationary and the integral equation ( 16 ) would be easily 
solved by the Fourier transform: kftjjt^) becomes with Fourier 

transform 



fjFWl 2 

T| +fjF(W )| 2 
o * 


(The tilda will denote Fourier transform throughout. ) 

; . ’ ' V* ,* 

This suggests a “quasi- stationary M solution form under (17): assume 
w ^ ere » Fourier transforming on the second argument 



k l = 


fjt ) | FW\ Z 

•V 1 

ri +f (t.) If(uj)| 

o * 1 




It may be verified by substitution into (16.1 using the assumed slow variation 
off with respect to F, that (18)is an approximate solution. 

It is to be expected, in view of the discussion in [l], that if dispersed 
transmitted pulses are used, the first step upon reception would be pulse 
compression and this data is then available for further processing. If the 
compressed pulse support is small relative to the distance over which g 
changes, we can then reasonably apply the results of this case. 

Realization - The quadratic functional Q(Z) given by Eq.15 can be 
written in a more special form now using Eq. 18 

i r • — r -m~ f * (t > l?(w)l * 

Q(Z) = ^ dt Z(t) J dtu e Z t { w) - .".T 


o T 


TL+ f # (t)|F(«))| 


^ 2 

(Z^,{t) s Z(t), tG T, and zero otherwise. ) Suppose | F(ttJ ) |- = 2 tt£^/ 0^, = P , 


u|< 0/2, and zero otherwise; then 

Q(Z) = J dt ZCt) (Z *h t ). 

o T 




T1 /P+f*(t) 

o + 


where h^ = 1, |u)j<n^/2, and 0 otherwise. It is likely the originally received 
data is already approximately so filtered and hence Z^^h^: if, e. g. , 

filtering matched to the transmitted pulse had already been done then 


^S-=-Z *h ; with this assumption 
. t t t 

Q(Z)=-i- [ dt|Z(t)l 2 . 
*0 J T 




T1 /P+fJt) 

o - p 


(15a) 


which is a weighted sum of the received instantaneous power after filtering 
matched to the transmitted pulse. 



By definition the .maximum likelihood estimator ( 11 MLE' 1 ) of the unknown 
parameter 0 is that 0 maximizing the log likelihood ratio which here has 
the form {-0 +Q(Z)} where [2, p. 177] presently we have 




ae. - -n 

i o TT 


II R z ( VV ae". k <VV dt i dt ; 


= 2^ J_ dt I dwff * (t > ^> + V • TT. k(t ’ “ > 


'o T 


The necessary conditions for the MLE's 0^ are found by setting the derivative 
S/S0. of the log likelihood ratio equal to zero: 




_ r 

Cn /P+fJt )] 2 

o *«• 



- [ri o /p+f 4 (t)]} . 


It is well known [4] that the correlation matrix of 
efficient estimates - those unbaised estimates of minimum possible error 
variance whose performance can be achieved asymptotically by MLE’s - is 


the inverse of the matrix C of elements 

c ij 3 - E {iSr A (z) } ' ’ 

where A(z) is the likelihood ratio and 0 = (6^ . . . , 6^} is the vector of 
unknown parameters. In the present problem (Gaussian signal in white 
Gaussian noise) it may be shown [2, p. 179 et. f. ] that 


C.. = 


*1 . 




dt. 


-L 

ae. 


R 


(t i*V 


90. k ( t 1 - t 2 ) - 




f * ( V 9 W and k(t rV = k i ( VVV : the change of variables 

( T St - t , u = t ) followed by a Fourier transform on the T variable gives the 
12 1 

representations iyt.ui) = f Jt) ? (« ) and as given by Eq. (5), 

yielding the form . .. 

c ij = ^r I dt I d “^ ) 

ij \ i ) 

where ^(uj) = |F(w)| 2 . Recalling only depends on 0, carrying out the 
indicated differentiations gives 

2 


1 

C.. = -^7 

ij 2 tt 


I dt J du, [ 


F<u>) 


^ 2 

t - n + f Jt)|F(w)| 

O ■ * 

To proceed assume that, reasonably, 

. i 2 


2 Bf j- 

i * 


af» 


a e. 


3 0 . 

3 


u)| < n F / 2 , 


. „ 7 f 2" V n F . 

. |F>)| 2 = { F F 

o , M>cy2 , 


( 20 ) 


where £ is the energy of the complex modulation (which is twice the energy 
F 


of the real, bandpass transmitted waveform). Then 




C.. = 


ij 2 tt 


i 


dt 


(Bf* /B0.)(Bf* /30.) 

t i t j 

C-T| /P + f*(t)] 2 
o 


tf9a) 


0 =2tt£ ]n . "We now examine the forms (20): they are somewhat involved 
F F 

though in specific instances numerical integration is straightforward. 

As an example of a useful approximation, note f^ results from the 
convolution of, essentially, the un variate density function of the gross height 
variation and the antenna pattern, as scaled into delay. A reasonable 
approximation in the denominator of (20)is arrived at by noting that, ignoring 
the sidelobes of the antenna pattern, we can take f^(t) as "pulse-like 11 , 
approximately of constant value - say g^f(0) (Note (2Q)does not depend 


on T ). 


fmt 

Hence (6) becomes setting B = fl / 2tt and s B ./(T} /£ + g f ) » 

JT JT Jc o o o 


C.. 


. &f* 

-^he- ■ it 


dt 


_ & r Sf *w 


; 


fcf. 


2tt j ae. 




30. 

J 


where fjiu) =g Kq, (w) cp [u> /(c/2)]exp (-iU) T ), and Cp is the univariate characte 
o 1 h o h 


function corresponding to p . 

. h 


-Consider first finding the covariance matrix of the efficient estimates of 

amplitude, delay, and sea state; assume the uni variant distribution of the 

2 

gross waveheight is normal with zero mean and variance then 

“ )Z /2 (c/2) 2 - i<UT 

= g 0 Kqj(»-') e 

~ - iu)T 0 . .. 

= g Q f (UJ ) e 


Define CT = a /{c / 2) and the vector of unknown parameters 0 ~ (T , Q,g ), 
fa o o 

Then we calculate: - • 


C,. =Ctg - M(l) , 

i-i . o 


C 12 = C 1 3 = ° ’ 

C 22 = • M(2) , 


C,, = -<Jg a M(l) , 

23 o 

C 33 = Qm(0) • 


1 f 2k ~ 2 

Here we have defined M(k) — J w / £ ( tu ) / dtu; of course M(0) = is the 
-energy of the function f and, setting M(l) -6 f , we may interpret Q f as the 



radius of gyration measure of the bandwidth of f^ (provided the centroid of 
|lj ^ is zero). Inverting this matrix we find the variance of the efficient 

j 

estimators: 


* 2 . 
Var(T ) = r— 

° CTg o Z M(l) 


* 1 1 
Var (o ) = j- * j— r 

1-p CIg o a M(2) 


Var(g*) = 



1 

Om(0) 


C’ 2 ‘la) 

(2'lb) 



121 c) 


where the correlation coeficient o of the efficient estimates of sea state and 
amplitude is 

C 2 2 ^ . 

p 2 - 23 _ M(l) ^ f21d) 

C 22 C 33 " J 

# . . * 
the efficient estimate T of delay is uncorrelated with the other two efficient 

estimates. 

■ The zero correlation between the delay and sea state efficient estimates holds 
more generally. Suppose = Cp^ 0 ^): fc k en correlation is zero ^ and 

only if the univariate distribution of the gross sea height is even. For then 
Bcp^(u)/Ba = ucp^u) and , 


0) 


U) 


\tzfz' \ :^ 7 ?> d " ■ 


12 


2tt 


~ ,2 


First note that, as q^ is real, | qj is the sum of the square of an even function 
and the square of an odd function and therefore even. Second, as p h is real, 



Cp is real and even if and only if p is even; more generally Re{cp cp 3 is 
n h h. h 

odd - and hence always contributes zero integral - and dtu{cp cp } is even 

h h 

and hence can make Cj2 non “ zer0, 

, ■ • - ■ We will be able to carry out the integrations 

required in Eq.(19a) above if we assume an exponential form for and cp^: 

let q^W) = Iq^O) | exp (- ^ u;| ) and cp fa (uu ) = exp (- a | U)| ). Then 

or. 


f * (t) = g o * .¥? 


[l + (±L T o_) ] 
* . J 


where ^ " t }>2 + ° * CT an d the parameter or = = K.|q(0) | is 

given explicity in terms of radar system parameters by Eq. 22 of [ l] : 

3 2 2 

or = AqQ/(4TT) R , where o =rr g /X . Also from [l], App. Ill, 
o o o o fr ’ 

2 , 

*j>j w X R^/ Ac. The elements fcyj given by(2d)may now be straight- 
forwardly, albeit somewhat laboriously, calculated by residue integration: 


33 


12 


°F 

r l *o a i* 

2tt 

\ ,2 

• V ■ 


V, 

■g a p 2 
h o 1 N 

2TT \ 

i 2 ' 

tt.* n 

o 


( "l e ) Z 

2rr 

\ TT ^ / 

C = 0 
13 

i 

2 2 

( 

g o“l P 


4> * 


JT 

2 


2 a(a+l) 3 


3 2 

a -a + 3a + 1 


a 3 (a+ l ) 3 


* * IT 


23 2rr 


2 3 2 
TT ,<j> n 


\ ; TT . 

)• ♦ ’T 


TT • a - 2a- 1 


3 2 * 

a (a+ 1) 


. 2 


where a = 1+g or. 0 /tt 6r\ ; recall 2 tt£ /0 . 

o 1 *- T . o F F 



#v 

Defining a ’= g a P 0 .^/tt < j> D = 2g or C / <}>D » we may write the 
o j. F^ o o 1 F o 


{c..} in the form 
1J 


C 11 = $ 


a a - 1 


a Z -l 


a(a + 1 ) 3 ^ a(a+ 1 


_ a- (a^ -l)(a 3 -a^ + 3a+ 1) ' 


22 4<j> 


a 3 (a+ l) 3 


O •> 

_ ; a_ ^ (a -l)(a -a + 3a + 1) 

4<J>’ 3, r,» 2 

Y a (a+ 1) 


_ a (j> a — 1 

33 ” 2 * 3“ 

4g o a 


C 1Z =C 13 = ° * 


and 


a (a 2 -l)(a 2 - 2a -I) _ a (a - l)(a 2 -2a -I) 


23 4g 


a 3 (a+ l) 2 


4g o a 3 (a+l) 


Inverting this matrix to find the covariance matrix of the efficient 
estimators, the diagonal elements are 


Var(T*) = -*■ , 

° a a-1 


(21e) 


. d> (a+ 1) 

Var (a ) = -* • 


- 1 


(21f) 


and 


* g c 
Var(g ) = r 
o ~ 


* 


a 3 -a^-t 3a + 1 
a-1 


(21g) 


As an application of these results, consider the choice of best pulse 

modulation bandwidth , which enters these error variances only via 
2 ^ 

a = 1+ a/Q^,: a can equivalently be chosen, then, to minimize one of these 



i 


three forms. It is easily found algebraically that a= 2 minimizes both 

' Var(C7 ) and Var(g ) and a = (3 + /l7)/4«L 8 minimizes Var(T ); further, a 

o o 

* 

choice of a= 2 increases Var(T ) by less than 2%. Thus the choice a = 2 is 
* o 

^ 2 2 
a very good one; then a/0 = (2TT f (0) C /t| )/2ttB = Q (0)/a is seen to be 

F * F o F s n 

(very nearly simultaneously) optimally 3. 

This choice of a is allowable, provided, roughly, the signal-to-noise 
ratio is large enough so that the quasi- stationary approximation remains 
valid: for this we must have 1/0 << 0. Since 0/0= 3 (for a= 2), 

X 

1/0= 3/a = (3/2)(r] /<*£• )0: thus, more precisely, we must have 

o F 


at 


< < 1 . 


As will be seen below, this condition is easily met in practice even for a 

satellite vehicle. — — - — -r — ' _ 

With a = 2, we find . 


Var <V = 9 


and 


* 27/^0 \ 2 

Var <° > = T ( rf" ) ( V a) • 


Var(g Q ) 


g. 


11 


(±h 


the correlation coefficient of the efficient estimates of sea state o and 
amplitude is, in magnitude, 1 //33, a rather low value in that Var(o ) and 

Var(g *) are increased by this correlation only by the factor (33/32). 

o 


(22a) 

(22b) 

(22c) 



Evaluating the approximate forms given by Eqs. 6 for this particular 


f we find 


VadT*) - 16 (^f- ) <V > 

F 

* 64 / \ .2 

Var(a ) « — (^-^r J (<l» 1 + °) 


(23a) 

(23b) 


Var(g 


o 


- PM - 

3 Ufi / 




2 3 

g o F 


(23c) 


and p 2 = 1/4; here <X = CX = tt£„ 4>/4T) or corresponding to unity SNR was used. 

opt -T o 

The parameter dependencies enter identically in the exact and approximate 
forms. Comparing the error standard deviations, the numerical factors, 
in order, are in the exact. calculation (3, 3. 7, 2. 4) and in the approximate 
calculation (4, 4. 6, 3. 3): these are in good agreement. 

A Specific Exact Solution. - Suppose that the reception time is large 
relative to the time over which the returned pulse is non-zero, a good 
approximation in practice: then in Eq. 4, the integral equation that determines 
the kernel function k(t^, t^) of the optimum processor, we may set T = (-»,»)• 
Taking the two-dimensional Fourier transform (denoted by a double tilda) of 
this equation we find 

*« 

T3 o k(u,v) + —J R s (u, -tt> )k (UJ , v)dw = R g (u, v) , -»<u, v < 


(24) 



From Eq. 8 of [l] we have (temporarily set T q = 0) 

Rg(u, -v) = p (u-v)F (u)F (v) 

where 

p (u-v) =K q (u- v) cp ) • 


Again assume 


2 rrC 


F(u) = 


0 


I u |< n F /2, 
|u|> o F n ; 


then Eq. 22 becomes 


y k^u, v) +. | pT (u-«) )k , v)do> 
■°F /2 


n 


p(u-v),u, v€ (- — 



(24a) 


(24a) 


where k^(u, v) = k (u, -v). 

The integral equation (22a) satisfied by k l is of a type well known in 
communication theory (see,, e. g. t 3] ). There is a straightforward method 
of solution when p is the (one-dimensional) Fourier transform of a rational 
function; thus any reasonably behaved q * p^ can be approximated by a 


rational function and the corresponding k found - in principle, though the 


calculations may be tedious. ■ _ 

To apply this method here assume that 

p(u) = a , -» < u < » ; 

2 

then p(t) = (a/rr<j>)r 1+ <t/ ) ] which roughly approximates a possible q * p^: 
it would be desirable perhaps to model the usual asymmetry but this requires 
more complication. The solution to Eq. 4b is now easily found by adapting 



the solution for an interval (0, 0^ found in, e, g. , [ l] , p. 388, Eq. 1. 46: 


l(u,v)= { 


K(u) K(-v) , -0 • /2 < u < v < 0„/2 , 
x x 

K(v) K(-u) , -0 /2 <v < u < n F /2 , 


( 25 ) 


where 


A P(n F /z+u) A -p(n /2+u) 

K(u) = ^2) [ (p+ <|))e + ( 3 - c|> )e ], 

2 a a 2 * 2 

2D e «{>orP^C ri 3 T(P + <j>)e -<P -<{,) e ' ] } , 


and |f s +2ptf<j>/ T l . Of course k itself is now readily found and further 
one can show (by an easy argument using the integral equation) that when 

T o* °* -IT (u+v) 

- ' ° -0/2 < u < v < 0 /2, 

F x 


(U - V) = { -iT (u+ v) 

e ° K(-v)K(-u), -0 F /2 < v < u < 0^2. 


Suppose P0 F >>1: then 

k(u,v) « • e- g|u+V ' ,- V 2<u,v< a n. 

v. o ? ( 9 + Pf 

A 2 2 2 2 4tt£ f“ 

Since p 0^ = 4» “p + “ * * ** * s sufficient for this case to 


o 


obtain if either 


4TT * — > > 1 - 

T F .2 


that is, the smeared returned pulse is long relative to the transmitted pulse 
resolution - or 


Tl. 


<j,a„ = 4 tt • 


2n£ F g _j>_ 

*n ‘ p -r-i 


>> i 



that is, the product of SNR and $/ 9 F is at least the order of unity. (In 

practical designs contemplated these conditions are met.) 

We know from [l] that the calculation of Q(z) is realized by the 

configuration shown in Figure 1 - because k has the some functional form 

as R which is k (aside from the multiplier T) in the threshold case 
b O 

which was discussed there. The final linear, time -invariant filter has a 
transfer function 


k(UJ ) « 


4>a 3 


ti o P(p+ <t,) 2 


• -Pi uul 

* e , - 05 < it) < ® : 


in particular its bandwidth is about 

r 1 « <— ) 1/2 ■ 

P K 2a$$ } 

The result agrees with the same results obtained under the quasi- stationary 
approximation. 

It would be of interest to use this precise result for calculation of 
performance limits and hence the best bandwidth setting for any SNR. 


Nonsymmetric Wave Height Distribution. - Although it is often 
reasonable to assume a normal univariate density function for the random 
gross wave height, it is known even in simple cases that the density function, 
is only approximately symmetric. It is clear that a nonzero mean can be 
completely ambigous with the geo id parameter: see Eq. (1).’ Such a 

nonzero mean may be associated with a large area current. We therefore 
consider a zero mean, nonsymmetric distribution that is a small perturbation 
from normal: the Gram-Charlier series f 3] when truncated gives a 



convenient representation:., thus assume 

-a t n 

V (t) = e (1 - i0 4a t) 


where 0 = E {h ]/3l will generally also be an unknown parameter. We 

therefore consider the efficient joint estimation of 0 = (T , CT, g , 0 

o o 4 

_ 3 3 

0 . = 0 , /(c/2) . It will be convenient to write cp, (t) = cp, (t) (1 -i0 A t ) where 
4 4a ha h 4 

2 2 

Cp (t) as before is the normal characteristic function (exp (-a t /2) and 
h h 


keep f and f^ unchanged. 


The matrix of c. ^ ' s is straightforwardly calculated to be the symmetric 


matrix 


C = <t g 2 


M(l)+0 M(4) 
4 


. 0 
0 

M(2) 


a 2 [M(2)+ 0 2 M(5)] 


— [M(l)+ e.M(4)] — , [l+e,M(3)l 

g 4 2 4 

° 


0 


- © 4 a M(4) 


M(3) 


M(3) 


Note that if the true (unknown) value of 0 , is zero then C becomes 

4 


M(l) 




a 2 M(2) 


g. 


M(l) 1/g 


M(2) 


0 


0 


M(3) 



Then the correlation matrix of the efficient estimates is 



_J 

M(l) 


1 


1 


1- P 


23 

2 


M(2)ct 


•g o P 2 “/M(l)a 


1 * P 


23 


1- ft 


2 


23 



1- P 


14 


1 

M(3) 


where 

. 2 , M(2) 2 
14 M(1)M(3) 


and 

2 _ M(I) 

P 23 M(0)M(2) 


The variances of the efficient estimates are, of course, the diagonal elements. 

In the more interesting case when the true value of 0^ is not zero we 

2 

may use the assumption that 0^ is small by ignoring the 0^ terms to ease the 
labor of calculating the inverse C * which may then be approximated as 
C 0.C l C,C ^ where 

o 4 o 1 o - . 




It is found that * G^ has zero diagonal elements and hence the variances 

of the efficient estimates are unchanged; that (1,4) and (2, 3) elements are 

zero and hence the correlation of the efficient estimates T and 0 , are 

o 4 

* * 

unchanged as with the efficient estimates q and g Q ; and that the (1, 2), (1, 3) 
and (2,4) elements are nonzero so that there is at least a weak correlation 
between all efficient estimates. 

Summarizing, 


* 1 

Var{T ) = j 

O j. 2 




Var( V= ~2 


<£g 


o 


1- P 


14 


1- P 


14 


1 


M(l) 

J 

M(3) 


and 


Var (a ) = 


<tg 


o 


Var (g Q ) 




1- P 


23 
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Example - Suppose that |q(U))| = ] q{0) | . exp(-c(> UJ ) so that 
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where = K| (0) | and v =[2(a+<J> )] . We then have M(0) = /2tt’v , 


M(l) = /2 t?v • v , M(2) = 3 ot*/Zn v * v 4 f and M{3) = 15 /2tt v * V 
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The cCrrelation coefficients are then p = 9/15 and p = 1/3 and 
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(1 - P ) = 15/ 6 and (1 - p ) = 3/2: these are the factors by which the 
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various efficient estimates variances are increased by their correlation. 



If the signal pulse bandwidth is chosen optimally we have seen 
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These results can also be written in the form (in part using the 
definition of V ) 
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Threshold case. - In the " threshold case" (in which the estimates 
can have arbitrarily good accuracy) we can readily find the structure of the 
best estimators and calculate the (C ramer-Rao) lower bound on the 
estimator's accuracy. From the best processor calculates the form 


Q(z) =-^; q(-t)* p h (- f t)*lF(-t)*Z(t)| Z | 
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the realization of this calculation by a 'taixed integrator" processor is 
discussed in [1]. .We find, by direct calculation, assuming geoid and sea 
state unknown, 
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Here we also assumed that the univariate sea height distribution is normal 
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[4], of zero mean and variance J Rg is the two-dimensional Fourier 
transform of the signal covariance function R . Thus fl] 
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where, then F is the signal pulse spectrum (Fourier transform) and cp^ is 


the univariate characteristic function of the random sea height. 



We note immediately that since the left side of (25b) is real, the right 

side must be zero: that is, the efficient estimates of geoid and sea state are 

uncorrelated! Note also that only the modulus If! of the signal pulse 

spectrum enters, then, into the accuracies. 

If we assume the signal pulse spectrum modulus form given by (20) 
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and | §( u ) I = lq (0) exp(- <|> u / 2)| , then we can calculate 
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and, finally. 
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Y(a) = [ e X dx 


Note that since the matrix of elements c.. is diagonal we have, as 

2 

aforementioned, the efficient estimates *r* and (a )* uncorrelated and their 

n 

error variances are just the reciprocals of (26a) and (26c) . These forms can 
now be used to make inferences about good system design just as in [ 1] . 

For example, note that the pulse modulation bandwidth 0^ enters forms 
(26a) and (26c) only as a = vOp in the arguments of andY^^# resp. 



Uncorrelated estimates. 


We have, in the course of the above 
calculations, showed that the efficient estimators of geoid and sea state are 
uncorrelated when the univariate probability density of sea heigh is Gaussian. 


We now show that it is the property of evenness of the density function 
that is sufficient for this property of the efficient estimates. (We are still 
assuming the threshold case, ) 

Recalculating (3) using (5) for an arbitrary characteristic function 
cp (u) = $(a u), the required derivative = u Cp' (o^u) =u Cp^'Cu). 

Then we may write 
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where G (t» ) = I F(u)) ! 2 and H(u) = | q(u) ) CP U ' (“ JJ ) ^ • Set 1 ec l" al to 

i times the double integral; then by Parseval's theorem we may write 
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X = i J ! G(t) | H(t)dt. 


In what follows denote the real part of a complex- valued function F by F R 
and the imaginary part by F^. 

We first note that since, in view of its definition, G is real, G is 

2 2 2 

even and G is odd and therefore JG ! = G + G is even . Similarly of 

I * R t i t 

2 

real implies | qj even. 

We also note that the realness of p implies cp is even and 

n n R 

CD is odd (cd s 0 if and only if p is even); hence cp' is odd and cp' is 

n I n I n n R Hj 

even and so fm CD 9 is odd and f cp Cp 'V is even ({ cp c 0 '1 s 0 if and only 
nn R n n I n n I 

if p is even). Thus H is odd and H is even, implying H is purely 
n K 1 



imaginary - thus checking that I is in fact real - and it will be odd if and 
only if p n is even. 

Thus finally I = - \\ G(t) | 2 H^tjdt will be zero if is odd which is 
true if and only if is even. 



.REFERENCES 


U. R. O. Harger, 'Radar altimeter optimization for geodesy over the 

sea, " IEEE Trans, on Aerosp. and Electr. Syst. , Vol. AES-8, 
November 1972, pp. 728-742. 

2. H. L. Van Trees, Detection, Estimation, and Modulation Theory. . 

Part III: Radar-Sonar Signal Processing and Gaussian 
Signals in Noise. New York: Wiley, 1971. 

nd 

3. C. W. Helstrom, Statistical Theory of Signal Detection, 2 Ed. 

New York: Pergamon, 1968. 

4. M. S. Longuet-Higgins, ’In the statistical distribution of the heights 

of sea waves,” J. Mar. Res., XI, No. 3, 1952, pp. 245-265. 

5. R. P. Dooley, F. E. Nathanson, L. W. Brooks, 'Study of pulse 

compression for high resolution satellite altimetry, " Final 
Report, NASA Contract NAS 6-2241, May 1, 1973. 



APPENDIX 


It is of interest to review the situation when the scattering mechanism 
is not random so that a signal known except for delay and spread is received 
in white, thermal noise. The well-known theory [31 for this model gives, 


e.g. 
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It is immediately clear that if cp is real - which is true if and only 

h 

if p is even - then =0: that is, the efficient estimates of delay and 

r n 12 

spread are uncorrelated. We have already observed that Refi cp cp ' t = 

h h 

f co Cp 'h is even and zero if and only if p is even. 
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We can also calculate readily 
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where we have assumed cp (u>) = exp(-a U) /2), a normal characteristic 
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function. 


If we assume, for ease of calculation, 
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The essential dependence of these forms on Q is shown in Fig. A.l 
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It is clear that 0 a small integer of multiple of 1 /a essentially 
F h. 

achieves the maximum value in C and C 0 insofar as they depend on p . 
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In this the results agree with the results for geoid and sea state estimation 


in the threshold case. 
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Figure 5. - 


Figure 6. - 


Geometry of radar altimeter. The coordinates are (x, y, z), 
where the x and z directions are the vehicle velocity vector 
v and the local vertical directions respectively. The antenna 
aperture is A and the nominal antenna pattern B has nominal 
beamwidth {3 ; F is an instantaneous radiated pulse position. 
The geoid G and sea surface are separated by h(x, y); the . 
geoid G has least range z = R q at (x = 0, y = 0). 


Realization of optimum processor when preprocessing SNR 

is small. A is a bandpass filter matched to the complex 

modulation F, B is an envelope - squared detector, and C and 

D are lowpass filters matched to q and p , respectively. 
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Weighting function w(t) = f Jt)[n /$ + fjt)] 1 for specific 
f* as a function of x= f o (°)7(r| 0 /pf* 

The weighting functions f v^(t)} appearing in likelihood equations. 
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Dependence of (/. (a) on a, where a = 1+a/Q ; (see Eqs. 21e, 
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Dependence of y (a) and y ( a) an a = vO_. 
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Fig.1. Geometry of radar altimeter. The coordinates are (x.y.z), 
where the x and z directions are the vehicle velocity vector v and the 
local vertical directions respectively. The antenna aperture is A and 
the nominal antenna pattern B has nominal beamwidth /? ; F is an 
| instantaneous radiated pulse position. The geoid G and sea surface 
are separated by h(x,y) ; the geoid G has least range z - /? 0 at 
! {* *= 0, jr ** 0). 
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Fig. 2. Realization of optimum processor when preprocessing 
SNR is small. / is a bandpass filter matched to the complex modu- 
lation F, B is an envelope -squared detector, and C and D are low- 
pass filters matched to <7 and p k , respectively. 
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